STUDIES  OF  THE  COEFFICIENT  OF  VARIATION  OF  THE  MAGNITUDE  OF 

EEG  SIGNALS 


W.  H.  Round1,  J.  W.  Sleigh2 

'Department  of  Physics  and  Electronic  Engineering,  University  of  Waikato,  Hamilton,  New  Zealand 
‘Department  of  Anesthetics,  Waikato  Hospital,  Hamilton,  New  Zealand 


Abstract  -  An  analysis  of  the  variation  in  magnitude  of 
EEG  signals  in  various  frequency  bands  of  anesthetized 
patients  and  normal  sleeping  volunteers  was  carried  out.  The 
coefficient  of  variation  (CoV),  i.e.  the  standard  deviation/mean, 
within  10  second  epochs  was  found  to  be  quite  constant 
throughout  the  whole  of  the  EEG  recordings  and  was  typically 
about  0.46.  This  was  found  to  be  the  case  for  both  the  patients 
and  the  volunteers.  Histograms  of  the  magnitudes  indicated 
that  the  magnitudes  are  distributed  as  fix)  =  p.ve"av  functions. 
However  a  CoV  of  0.46  is  consistent  with  fix)  =  p.ve"ar 
functions.  The  non-stationary  nature  of  the  EEG  is  such  that  it 
is  likely  that  while  over  short  periods  the  EEG  magnitudes  are 
distributed  as  fix)  =  p.ve"ar  functions,  variations  of  a  over  time 
mean  that  in  the  long  term  the  EEG  magnitudes  are 
distributed  as  fix)  =  p.ve"ar  functions. 
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I.  Introduction 

The  EEG  is  has  long  been  used  as  a  diagnostic  tool  to 
analyze  brain  activity.  Various  models  have  been  suggested 
to  describe  the  way  that  the  brain  produces  the  EEG  [e.g., 
1].  The  EEG  has  a  chaotic  nature  and  is  produced  by 
sources  that  are  quite  non-stationary,  but  such  systems  are 
often  studied  over  sufficiently  short  time  intervals  that  they 
may  be  considered  stationary  [2,3].  Here  we  analyze  the 
magnitude  of  the  EEG  in  terms  of  its  covariance  in  various 
frequency  bands. 

II.  Methodology 

EEG  signals  were  recorded  for  107  patients  undergoing 
general  anesthesia  for  surgery  and  6  normal  volunteers 
during  normal  sleep  without  anesthesia.  The  data  was 
acquired  with  an  Aspect  A 1000  EEG  monitor.  The  EEG 
was  sampled  at  256  samples  per  second.  No  filtering  was 
applied  to  the  data  during  recording. 

Spectrograms  were  calculated  using  the  Fast  Fourier 
Transform  [4]  for  each  of  the  EEG  signals  using  a  256  point 
Hamming  window  with  128  point  overlap  of  successive 
windows.  This  resulted  in  the  spectrogram  being  calculated 
at  0.5  second  intervals  and  in  1  Hz  bands.  No  pre-filtering  or 
artifact  removal  was  applied  to  the  data.  The  spectrogram’s 
magnitude  at  a  chosen  frequency  was  divided  up  into  ten- 
second  epochs  containing  21  samples  of  the  magnitude.  For 


each  of  these  epochs  the  mean  and  standard  deviation  of  the 
magnitude  samples  was  calculated.  Also,  the  coefficient  of 
variation  (CoV)  (i.e.  the  standard  deviation/mean)  was 
calculated. 


III.  Magnitude  variation  and  anaesthetised  patients 

Graphs  of  these  parameters  for  the  20  Hz  frequency 
band  of  a  patient  undergoing  general  anesthesia  are  shown 
in  Fig.  1.  Also  shown  in  Fig.  1  is  the  associated  EEG  and 
the  spectrogram.  In  the  spectrogram  the  darker  regions 
represent  the  lowest  activity,  and  the  lighter  regions  the 
highest. 
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Fig.  1 .  The  EEG  and  graphs  of  associated  derived  parameters  for  a  patient 
undergoing  general  anesthesia.  The  top  two  graphs  are  the  EEG  and  its 
associated  spectrogram.  The  next  three  graphs  are  of  (from  top  to  bottom) 
the  mean  magnitude,  the  standard  deviation  and  the  coefficient  of  variation 
of  ten  second  epochs  for  the  20  Hz  band.  The  bottom  graph  is  the  spectral 
entropy  of  the  EEG  signal. 

The  vertical  dashed  line  represents  the  time  at  which 
the  anesthetist  judged  that  the  patient  was  unconscious.  The 
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data  collection  was  terminated  before  the  patient  woke  up 
again.  In  the  spectrogram  it  is  seen  that  the  patient’s  EEG 
activity  was  initially  high  at  all  during  the  early  stages  of  the 
data  collection  when  the  patient  was  awake.  A  collapsing  of 
the  activity  at  the  higher  frequencies  as  the  patient  went  to 
sleep  is  also  seen.  Also  seen  in  the  EEG  are  large  amplitude 
signals  at  approximately  3.5  min,  4.5  min  and  13  min.  They 
are  probably  EMG  artifacts,  and  they  can  also  be  detected  in 
the  spectrogram. 

Also  shown  is  a  graph  of  the  spectral  entropy  to  act  as 
an  indicator  of  patient  consciousness.  It  is  noted  that  the 
spectral  entropy  decreased  considerably  some  20  seconds 
before  the  time  that  the  anesthetist  judged  the  patient  lost 
consciousness,  thus  loss  of  consciousness  may  have 
occurred  somewhat  than  was  realized. 

In  the  graph  of  the  mean  magnitude,  the  decrease  in 
magnitude  as  the  patient  went  to  sleep  can  be  clearly  seen, 
and  at  this  frequency  the  magnitude  increased  somewhat 
from  about  9  minutes  to  13  minutes  which  can  also  be  seen 
in  the  spectrogram. 

The  similarity  between  the  mean  magnitude  graph  and 
the  standard  deviation  graph  should  also  be  noted.  They  are 
similar  in  form  with  the  peaks  in  the  magnitude  graph  also 
being  seen  in  the  standard  deviation  graph.  It  should  be 
noted  that  large  amplitude  artifacts  in  the  EEG  correspond 
to  the  larger  peaks  in  the  covariance  graph.  This  is  probably 
to  be  expected  as  the  nature  of  the  standard  deviation  is  such 
that  it  will  emphasize  contributions  to  it  from  the  values  that 
differ  most  from  the  mean,  and  larger  differences  are 
expected  to  occur  with  the  large  amplitude  artifacts. 

Generally,  though,  the  standard  deviations  values  are 
just  under  0.5  of  the  corresponding  mean  magnitude  value. 
This  is  true  for  the  entire  duration  of  covariance  graph,  even 
during  the  initial  part  when  the  EEG  activity  is  higher. 

Consider  now  Fig.  2.  In  this  figure,  covariance  graphs 
have  been  calculated  for  the  EEG  signal  in  Fig.  1,  but  now 
for  5,  15,  25,  35,  45  and  55  Hz  frequency  bands.  Often  the 
covariance  graph  peaks  in  any  one  frequency  band  are  also 
seen  at  other  frequency  bands.  Also  the  covariance  tends  to 
generally  be  just  less  than  0.5  (the  modal  value  is  0.46)  for 
all  frequencies.  This  suggests  that  if  histograms  of  the 
magnitude  values  in  different  frequency  bands  are 
constructed,  then  they  will  possibly  be  of  a  similar  form  as 
the  ratios  of  the  mean  and  the  standard  deviation  of  each  are 
similar. 

To  find  the  possible  form  of  the  histogram,  histograms 
were  constructed  of  the  magnitudes  for  each  frequency  band 
over  the  entire  time  of  the  EEG  recording.  These  are 
presented  in  Fig.  3  for  the  EEG  waveform  in  Figs.  1  and  2. 


Fig.  2.  Graphs  of  the  covariance  of  the  EEG  signal  in  Fig  1.  in  different 
frequency  bands. 


Fig.  3.  Histograms  of  the  magnitudes  of  the  EEG  waveform  in  Fig.  1  in 
various  frequency  bands. 

The  histograms  do  indeed  have  a  similar  form.  The 
shapes  of  the  histograms  are  such  that  they  can  be  described 
by  functions  of  the  form 

./(x)  =  p.vne“m 

where  x  is  the  magnitude.  When  x  is  small,  the  xa  term 
dominates  and  influences  the  initial  rise  of  the  histogram. 
For  larger  x,  the  e"a'  term  dominates  and  influences  the 
later  decrease  in  the  histogram,  a  also  influences  the  x  value 
at  which  the  function  reaches  its  maximum,  p  determines 
the  area  of  the  histogram. 


Noting  that  for  n  >  1  thaty(x)  is  concave  towards  the 
right  for  low  values  of  x  whereas  for  n  =  1  the  function  is 
convex  similar  to  the  histograms  in  Fig.  3,  it  is  reasonable  to 
expect  that  n  =  1 . 

For  n  =  1,  and  for  low  integer  values  of  m,  the 
coefficients  of  variance  of  distributions  described  by  such 
functions  are 

m  =  1,  CoV  =  0.7071 

m  =  2,  CoV  =  0.5227 

m  =  3,  CoV  =  0.4574 

For  n=l,  then  the  covariance  values  indicate  that  for  the 
typically  observed  covariance  value  of  0.46  that  m=3.  Thus 
the  distributions  should  have  histograms  of  the  form 

fix)  =  fke  “3 . 

However  coefficients  of  variation  for  the  distributions 
in  Fig.  3  are  rather  higher  than  0.46,  and  are,  on  average, 
0.62.  This  indicates  that 

1  <  m  <  2 

rather  than  3 . 

Choosing  m=2  and  fitting  functions  of  the  form 
fix)  =  xe"“'~  to  the  histograms  gives  curves  as  in  Fig.  4. 


Fig  4.  Dotted  curves  are  corresponding  graphs  in  Fig.  3.,  solid  curves  are 
fitted  curves  of  the  form  (3x  exp(-ocx2). 

It  is  seen  that  the  curves  are  generally  well  described  by 
fix)  =  pxe_0L'~  curves  except  in  the  lower  frequency  bands. 
The  differences  are  because  the  EEG  is  not  that  of  a 
stationary  system  and  there  are  longer-term  variations  in  the 
magnitude  at  any  frequency.  Hence  the  histograms  in  Fig  4 
are  the  summations  of  many  shorter-term  histograms. 
Shorter  term  histograms,  if  they  could  be  formed  with  just 
21  samples  taken  over  10  seconds  as  was  done  earlier, 


would  reasonably  be  expected  to  be  of  the  form^x)  = 
(3xe"a'  and  have,  on  average,  coefficients  of  variation  of 
0.4574. 

At  low  frequencies  where  there  is  a  lot  of  variation  in 
the  EEG  magnitude  (as  can  be  seen  in  the  spectrogram),  the 
longer  term  histograms  exhibit  long  tails  from  the  histogram 
summations.  At  higher  frequencies  where  there  is  not  as 
much  variation  in  the  spectrogram  and  the  system  is  more 
like  a  stationary  system,  the  tail  is  not  as  pronounced. 
However  it  should  be  noted  that  the  higher  frequency 
bands’  histograms  exhibited  two  peaks:  a  small  peak  at  large 
magnitudes  due  to  the  initial  high  activity  before  the  patient 
went  to  sleep,  and  a  large  peak  at  small  magnitudes  when 
the  patient  was  asleep.  The  smaller  peak  was  excluded  in 
fitting  the  functions  to  the  histograms  and  calculating  the 
coefficients  of  variation. 

How  well  a  summation  of  pxe"“'J  functions  can 
approximate  a  (3xe"ar  function  is  demonstrated  in  Fig.  5. 
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Fig.  5.  The  unbroken  curve  is  generated  from  a  summation  a  gaussian 
distribution  of  functions  of  the  form  (3x  exp(-ax3)  chosen  so  that  it  is 
similar  to  the  55  Flz  band  curve  in  Fig  4.  The  dotted  curve  is  a  fitted  curve 
of  the  form  px  exp(-ax2),  while  the  dot-dash  curve  is  a  fitted  curve  of  the 
form  [fv  exp(-ax). 

In  Fig.  5,  the  solid  curve  is  generated  with  from  a  sum 
of  many  pxe"a'  functions  with  their  peaks  evenly  distributed 
along  the  x-axis.  The  position  of  a  peak  is  determined  by  the 
a  value  for  each  function.  Their  areas  given  by  a  gaussian 
curve  centred  at  x  =  18  and  with  a  standard  deviation  of  8. 
These  values,  along  with  sum  of  their  areas,  are  chosen  so 
that  the  solid  curve  is  similar  to  the  curve  fitted  to  the  55  Hz 
band  data  in  Fig  4.  The  dotted  and  dot-dash  curves  are 
pxe'a'~  and  pxe"“'  functions  respectively  fitted  to  the 
unbroken  curve  so  that  their  peaks  coincide.  Note  that  the 
Pxe"ar  curve  is  close  fit  while  the  pxe"“'  curve  isn’t.  If 
attempts  are  made  to  fit  the  pxe"a'  curve’s  tail  by  altering  a 
so  that  the  rate  of  fall  off  is  similar  to  the  unbroken  curve’s, 


the  curve  moves  considerably  to  the  right  and  the  fit  is 
worse. 

Over  any  time  interval,  the  distribution  of  (Ire"'1'  curves 
will  not  necessarily  be  just  a  simple  gaussian  distribution.  If 
the  EEG  does  not  fluctuate  in  magnitude  to  any  great  extent, 
then  it  is  not  unreasonable  for  it  to  have  such  a  distribution. 
However,  if  the  magnitude  does  fluctuate  significantly  or  is 
contaminated  by  artifact,  then  the  distribution  is  unlikely  to 
be  gaussian  or  of  any  other  simple  function,  but  may  have 
multiple  peaks.  In  such  a  case  it  is  unlikely  to  approximate  a 
p.re"“'  curve.  At  times  when  larger  magnitude  signals  from 
artifacts  combine  with  other  shorter  magnitude  signals,  the 
standard  deviation  of  the  magnitudes  will  increase  more 
than  the  mean  of  the  magnitudes,  and  thus  the  covariance 
will  rise.  Thus  the  covariance  can  be  seen  as  an  indicator  of 
stationarity  of  the  signal.  The  more  it  rises,  the  more  the 
signal  is  fluctuating. 

IV.  Magnitude  variation  in  normal  sleep 

The  coefficients  of  variation  for  10  second  epochs  were 
calculated  for  the  six  volunteers.  Graphs  of  these  are  shown 
in  Fig.  6  for  one  of  the  volunteers,  along  with  the  spectral 
entropy.  Note  that  during  all  sleep  phases  the  covariance  is 
again  approximately  0.46  (modal  values  are  0.435,  0.456, 
0.453,  0.453,  0.476  and  .459  respectively  for  the  55,  45,  35, 
25,  15,  and  5  Hz  bands)  at  all  frequencies  (although  it  does 
go  significantly  higher  than  this  at  certain  times).  This 
suggests  that  the  histograms  of  the  coefficients  of  variation 
at  the  different  frequencies  have  the  shape  of  (3xe"a' 
functions  for  persons  in  normal  sleep. 

V.  Discussion 

The  data  and  analysis  reported  here  for  two  EEG 
recordings  is  typical  of  all  those  studied.  It  seems  quite 
surprising  that  the  CoV  of  the  EEG  is  so  similar  throughout 
all  of  the  EEG  recordings,  and  also  that  it  is  so  similar  at  all 
frequencies.  For  the  patients  undergoing  general  anesthesia 
and  for  the  patients  in  normal  sleep,  in  the  short  term  it 
tends  to  be  approximately  0.46  both  when  they  are  fully 
asleep,  when  they  are  quite  awake  and  at  various  stages  in 
between.  This  suggests  that  if  the  magnitude  in  any 
frequency  band  is  sampled,  the  samples  are  likely  have  a 
distribution  which  is  similar  in  all  frequency  bands 
regardless  of  the  state  of  consciousness  of  the  subject.  We 
have  suggested  here  a  possible  form  for  the  distribution,  but 
note  that  over  long  periods  of  time  the  signal  changes  as  the 
system  that  generates  the  EEG  is  non-stationary.  Thus  the 
distribution  will  alter  and  the  function  will  change,  but  it  can 
still  be  described  by  a  function  of  the  same  basic  form. 
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Fig  6.  The  EEG  and  graphs  of  the  coefficient  of  variation  of  ten  second 
epochs  for  the  various  frequency  bands  for  a  volunteer  in  normal  sleep.  The 
bottom  graph  is  the  spectral  entropy  of  the  EEG  signal. 
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